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Abstract 

Some recent results for a three-dimensional, semi-classical, tight- 
binding model for DNA show that there are two types of polarons, 
named radial and twist polarons, that can transport charge along the 
DNA molecule. However, the existence of two types of base pairs in 
real DNA, makes it crucial to find out if charge transport also exist 
in DNA chains with different base pairs. In this paper we address 
this problem in its simple case, an homogeneous chain except for a 
single different base pair, what we call a base-pair inhomogeneity, 
and its effect on charge transport. Radial polarons experience either 
reflection or trapping. However, twist polarons are good candidates 
for charge transport along real DNA. This transport is also very robust 
with respect to weak parametric and diagonal disorder. 
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1 Introduction 



Charge transport in biomolecules is of relevant in many biological functions. 
Moveover, electronic transport plays a well known and fundamental role in 
many of biological processes, such as nervous conduction, photosynthesis, cel- 
lular respiration and redox reactions. Recently, DNA-mediated charge trans- 
port has focused much attention with experiments that have been conducted 
to determine its conductive properties, and theoretical studies have addressed 
the possibility of charge transport and its efficiency p. The importance of 
charge transport through DNA is significant because some mutations in liv- 
ing systems and radical migrations are critical issues in carcinogenesis studies 
and may yield insight into damage prevention or repair processes |21 13 El • 
Moreover, charge transport has been demonstrated to proceed within HeLa 
cell nuclei pi as well as in the nucleosome core particles ^ , and can provide 
a practical method of genetic screening for known gene sequences and an al- 
ternative method to hybridization based arrays 0. Material scientists have 
thought that DNA has a fundamental physical interest for the development 
of DNA-based molecular technologies, as it possesses ideal structural and 
molecular recognition properties for use in self assembling nanodevices with 
a definite molecular architecture |H|. 

Photochemical electron transfer experiments have been carried out with 
femtosecond resolution, and have shown, with the precise control of the DNA 
sequences, a fast time scale (< 5 ps) process, attributed to charge mobility 
along an unperturbed double strand, and a slower process (< 75 ps), that 
they believe refiects charge hopping between perturbed domains along the 
DNA strand 0. 

Experimental results have explained that DNA acts as a linear chain with 
overlapping vr orbitals located at the stacked base pairs and its conductivity 
is very sensitive to disruption, caused by base-pair mistakes or interactions 
with proteins, in the base-pair stack. Depending on the base pair sequences 
the transfer can be slowed or inhibited [TU] . 

Moreover, the robust, malleable one-dimensional structure of DNA can 
be used to design electronic devices based on biomaterials [THIT^ITnilT^ IT^]. 

The structure of the bent double helix A-DNA can be modelled by a net- 
work of oscillators taking into account deformations of the hydrogen bonds 
within a base pair and twist motions between consecutive base pairs. The 
three-dimensional semi-classical tight-binding model for DNA that was first 
proposed in Ref. ^21 ^] niakes also the assumption that the electron mo- 
tion is predominantly infiuenced by vibrational modes of the double helix. 
The nonlinear interaction between the electron and the vibrational modes is 
responsible of the formation of polarons or electron-vibron breathers. These 
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are localized excitation patterns that can be static, producing charge local- 
ization, or mobile, bringing about charge transport. It has been considered 
some other models with inhomogeneities [TSj . 

A variant of this model has been proposed in Ref. ^Hl, which lifts the 
restriction of the previous one, that the perturbation of the spatial vari- 
ables were small enough to perform a linear approximation in the dynamical 
equations. Two types of polarons appear depending on the parameter a, 
that describes the coupling between the transfer integral and the distance 
between nucleotides, and the form of activation: radial polarons and twist 
polarons. For radial polarons, the deformations associated by the electronic 
charge affect mainly to the radial variables, and for twist polarons, to the 
angular ones. In the last case, the moving polaron is slower than in the previ- 
ous one, but it is more robust with respect to the introduction of parametric 
disorder. See reference above for details and Ref. jSOI for twist polarons in 
other model. 

In this paper we consider the movement of polarons along DNA chains 
made out of different base pairs. The most simple one consists of homo- 
geneous DNA chains except for a single different base pair, which we call 
base-pair inhomogeneity. 

We have found that when charge transport is mediated by means of a twist 
polaron, this moving carrier can transport charge along a base-pair inhomo- 
geneity in a very efficient way. Moreover, charge transport persists under 
the introduction of a high degree of parametric disorder, that represents the 
inherent disorder in the surrounding medium and the inhomogeneous distri- 
bution of counterions along the DNA duplex j2I] . On the other hand, radial 
polarons are either reflected or trapped by the base-pair inhomogeneity, mak- 
ing them poor candidates for charge transport along heterogeneous DNA at 
least in the framework of the tight-binding models. 

2 DNA model 
2.1 Description 

We consider a variant ^H] of the model Hamiltonian for charge transport 
along DNA that was first proposed in JE]. This variant consists of con- 
sidering the full dynamical equations of the system instead of performing 
linear approximations based on the assumption that the deformations of the 
molecule were small. 

These models were designed to implement the basic characteristics of the 
DNA double helix structure needed for a proper description of the charge 
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transport dynamics. The starting point is the twist- opening model [221 12HI 
121] that has taken into account the hehcoidal structure of DNA and the 
torsional deformations induced by the opening of the base pairs. The bases 
are considered as single nondeformable objects. The helicoidal structure 
of DNA is described in a cylindrical reference system and the n-th base 
pair has two degrees of freedom, namely (r„,0„), where r„ represents the 
radial displacement of the base pair from the equilibrium value Rq, and (pn 
represents the angular deviation from equilibrium angles with respect to a 
fixed external reference frame. As we are interested in base pair vibrations 
and not in acoustic motions, the center of mass of each base pair is fixed, 
i.e., the two bases in a base pair are constrained to move symmetrically 
with respect to the molecule axis. Moreover, the distances between two 
neighboring base pair planes will be treated as fixed, because in the axial 
direction DNA is less deformable than within the base pair planes [22] • 
A sketch of the helicoidal structure of the DNA model is shown in Fig. ^ 

The transversal displacements of the base pairs are deformations of the 
H-bonds and the angular twist and the radial vibrational motion evolve in- 
dependently on two different time scales, then they can be considered as 
decoupled degrees of freedom in the harmonic approximation of the normal 
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mode vibrations ^3] . 



2.2 Hamiltonian 

The Hamiltonian for the electron transport along a strand of DNA is given by 
H = Hei + Hrad + Htwist, where Hei corresponds to the part related with the 
particle charge transport over the base pairs, if^ad describes the dynamics of 
the H-bond vibrations and Hfwist is the part corresponding to the dynamics 
of the relative twist angle between two consecutive base pairs. This electronic 
part is described by a tight-binding system of the form: 

Hei = ^En\n){n\ - 14-i,n|n - - Vn+i,n\n + , (1) 



where In) represents a localized state of the charge carrier at the n^^ base pair. 
The quantities {^,ri-i} are the nearest-neighbor transfer integrals along base 
pairs, and {En} are the energy on-site matrix elements. A general electronic 
state is given by = J2n'^nit)\n) , where c„(t) is the probability amplitude 
of finding the charged particle in the state \n). The time evolution of the 
{cnit)} is obtained from the Schrodinger equation ih{d'^ /dt) = Heil"^) 

The nucleotides are large molecules and they move much more slowly than 
a charged particle, then the lattice oscillators can be described classically and 
Hrad and Ht^ist are, de facto, classical Hamiltonians. For homogeneous chains 
they are given by (omitting the hat over them): 



H 



rad 



H 



twist 



E 

n 

E 



2J^^"^ 2 



(2) 
(3) 



where pJi ^ind pf^ are the conjugate momenta of the radial and angular co- 
ordinates, respectively. In these expressions M is the mass of each base pair 
(M = 2m, being m an average estimation of the nucleotide mass), J = M 
is the moment of inertia of each base pair and Qr,n is the linear radial fre- 
quency which is proportional to the strength of the hydrogen bonds. Indeed, 
for an homogeneous chain all these frequencies are equals, that is, „ = bnflr 
with bn = I Vn. 

In this paper we are interested in the study of charge transport along 
DNA strands where all the base pairs are of the same type, say A-T (or 
C-G), except only one of them which is of a different type, say C-G (or A- 
T). We can suppose that the ratio between the elastic constants of bonds 
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in a C-G base pair and an A-T base pair is 3/2 because the first involves 
three hydrogen bonds and the second two of them. Then, as in Ref. [211 
we take 6„ = 0.8 for an A-T base pair, and 6„ = 1.2 for a C-G base pair. 

is the hnear twist frequency, and we represent by 9n,n-i = (0n — <Pn-i), 
the deviation of the relative angle between two adjacent base pairs from its 
equilibrium value Oq. 

In general, the ionization potential of different nucleotides differs an 
amount of 0.2-1.0 eV j22j. In our model, this implies different values of 
the on site energies for each base pairs. In this work, we have focused in 
geometrical effects due to the stretching of the chain over the charge, and we 
have considered E'^ independent on the type of base. 

The electronic part of the Hamiltonian, H^i, has a dependence on the 
structure variables r„ and 0„ through the dependence of the matrix elements 
En and Vn,n-i on them. The energy on-site matrix elements are given by j2Hj 
En = E^ + kvn, expressing the modulation of the on-site electronic energies 
{E^} by the radial deformations of the base pairs. The transfer matrix 
elements Ki,n-i; which are responsible for the transport of the electron along 
the stacked base pairs, are assumed to depend on the distances between two 
consecutive bases along a strand dn,n-i as Vn^n-i = Vo(l~«'^n,n-i)) where a is 
a parameter that describes the influence of the distance between nucleotides, 
and the later is determined by 

dn,n-l = [a^ + {Ro + Tn)^ + {Ro + Tn-l)^ - 

2{Ro + r„)(i?o + cos(^o + On,n-i)Y^'' - (4) 

with /o = (a^ +4/^2sin2(^o/2))^/^ where a is the vertical distance between 
two consecutive base pairs. 

In this paper, as in Ref. |1£,, we have not used the expansion of this ex- 
pression up to first order around the equilibrium positions, as was performed 
in Refs. [T7] . This allow us to consider parameters that allow larger 
deformations in the angular variables. 

Realistic parameters for the DNA are given in Refs. [23 1201 • We have 
considered: a = 3AA, m = 300 amu, Rq = 10A,fi,. = 8 x 10^^ s"\ Q^p = 
9 X 10^^ s~^. Ab initio calculations find that between adjacent nucleotides, 
the transfer integral is in order of 0.1-0.4 eV [SO] . Although this transfer 
integral is different for each pair of different nucleotides, we will consider the 
same value for all neighboring cases Vq = 0.1 eV, a supposition widely used 
and can be valid in order to reproduce ab initio results and experiments [SI] . 

We scale the time according to t —>■ Qrt, and we introduce the dimen- 
sionless quantities: f„ = r„ {MQf./VoY^'^ , K = K/ {MVLlVof^ , En = K/K), 

n = n^/nr, v = Vo/{jnl), a = a {Vo/M^ilf/\ Ro = Ro {Mnl/v,fi\ 
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2.3 Dynamical equations 



Using the expectation value for the electronic contribution to the Hamilto- 
nian, the new classical Hamiltonian H = {(f)\H\(f)) /Vq is given in the scaled 
variables (omitting the tildes) by 



n 

+ krn)\cn\^ - (1 - a(i„,„-i)(c*c„_i + c„c*_i) [. (5) 



The Schrodinger equation and the Hamiltonian equations lead to the scaled 
dynamical equations of the system: 

irc„ = {El + krn)cn 

- {I - a dn+l,n) Cn+l - {I - a dnn-l) Cn-l, (6) 



— a 



i>n+l 

aV 



,(7) 



-fi^ (20„ - (f)n-l 



^-7 (c„C„_i + C„C„_ J + — (c„+iC„ + Cn+lCj 18) 



where the quantity r = nQr/Vo measures the time scale separation between 
the fast electron motion and the slow bond vibrations. In the ordered case, 
with En = Eq, Wn, with the limit of a = the set of coupled equations 
represents the Holstein system, widely used in studies of polaron dynamics 
in one-dimensional lattices. Also, for a = k = 0, and random E^, the 
Anderson model is obtained. 

The scaled parameters take the values r = 0.053, = 0.013, V = 
2.5 X 10~^, Rq = 63.1 and Iq = 44.5. We fix the value k = 1 and consider the 
parameter a as adjustable. 



3 Charge transport mediated by mobile po- 
larons 

3.1 Mobile polarons 

The principal interest of this paper consists in the study of charge transport 
along the double strand by moving polarons in presence of a base-pair inho- 
mogeneity. We need first to obtain localized stationary solutions of Eqs. (jEllHI)- 
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In Appendix A, we present the procedure that we have followed for obtaining 
these stationary solutions applied to different cases: homogeneous chains, and 
inhomogeneous chains. The polaron motion can be activated under certain 
conditions. A systematic method to do this, known as the pinning mode- 
method [22j|, consists of perturbing the (zero) velocities of the ground state 
with localized, spatially antisymmetric modes obtained in the vicinity of a 
bifurcation. This method leads to moving entities with very low radiation, 
but has the inconvenience of being applicable only in the neighborhood of 
certain values of the parameters. An alternative is the discrete gradient 
method [SBl, perturbing the (zero) velocities of the stationary state {r„(0)}, 
{0n(O)} in a direction parallel to the vectors (Vr)„ = (r„+i — r„_i) and/or 
(V(/>)n = (0n+i — 4>n-i)- Although this method does not guarantee mobility, 
it nevertheless proves to be successful in a wide parameter range. We will 
denote the energy of the perturbation as AK = where A is the modulus 
of the vector used to perturb the system (parallel to Vr„ and/or V0ri)- This 
energy gives us the difference of the energy between the moving polaron and 
the static one, and it is usually called the activation energy. 

3.2 Homogeneous chains 

In this section we recall the basics results obtained in Ref. ^|, see this ref- 
erence for details. For homogeneous chains in absence of diagonal disorder, 

= E^,'^n, and if the parameter a is small enough, it is not possible to 
move the polarons by perturbing the angular variables. Mobility can be ac- 
complished only by perturbing through the radial variables. However, for 
larger values of the parameter a (a > 0.01), the polarons can only be moved 
by perturbing the angular variables, i.e., perturbations of the radial variables 
cannot activate mobility. Nevertheless, for intermediate values of the param- 
eter a (0.005 < a < 0.01), the polarons can become mobile by perturbing 
any set of variables, the radial or the angular ones. The movement is rather 
different, when the polaron propagation is activated by means of only radial 
perturbations, its characteristics are similar to the radial movability regi- 
men, and likewise when it is activated by angular perturbations. A detailed 
analysis of these different regimes can be found in the reference above. In 
general, radial movability requires less energy and has higher velocity than 
the angular one. Also, the limits of these regimes are not exact, depending 
on the kinetic energy of the perturbation. 

If an amount of disorder in the on-site energies is introduced, with 
random values l-ElJl < AE, we find that moving polarons exist below a criti- 
cal value AEcrif Beyond this value, polarons cannot be moved. In general, 
as shown in Fig. |2l corresponding to the mixed regime where polaron can be 



8 



3,5rT 




0,01 0.02 am OM 0.05 0,06 O.O? 0,08 

i K (eV) 



Figure 2: Velocity of the polaron as a function of the energy amount AK 
corresponding to the mixed regime when the polaron is activated by means 
of angular perturbations (a = 0.01). Circles and solid hues represent the 
disordered case with AE = 0.05 and triangles and dashed line the ordered 
case. 

activated (in absence of disorder), by means of both angular and radial per- 
turbations, the mobility induced by angular activation is more robust with 
respect to parametric disorder. If disorder is high enough, radial perturba- 
tions that could move polaron destroys it. Thus, it only can activated by 
means of angular perturbations. In this case, the movement is similar to 
the ordered case, the polaron has lower velocity, and the activation energy is 
higher than in the radial movability regime. 

As shown in Fig. |21 it can be appreciated that for a polaron in the mixed 
regime, and if AE is high enough, it is impossible to move it with radial 
perturbations, but only with angular perturbations. In this situation, the 
movement is very similar to the ordered case. 

3.3 DNA chains with a base-pair inhomogeneity 

The studies of polarons in homogeneous chains are only applicable to syn- 
thetic DNA made out of a single type of base pair. In real DNA, the two 
different base pairs, A-T and G-C, combine in different ways constituting 
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Figure 3: Soft inhomogeneity and radial polarons. Trapping phenomenon 
due to the interaction between a moving radial polaron in a C-G chain with 
an A-T base pair 

the genetic code. Thus, as a first step in this direction, we consider a ho- 
mogeneous DNA-chain except for a single base pair of a different type. Sys- 
tems of this type are easily synthesized and the numerical and theoretical 
results could be compared with the experimental ones. The outcome would 
help to determine the actual DNA parameters and, therefore, the type of 
polarons that can be expected. There are two different systems: a) a hard- 
inhomog entity, that is, an A-T DNA chain with a C-G inhomogeneity; b) a 
soft inhomogeneity, i.e., a C-G DNA chain, with an A-T inhomogeneity. 

We activate the motion of a stationary polaron, using the discrete gradient 
method, at a site about a hundred sites far away form the location of the 
inhomogeneity. Different types of moving polarons can be obtained : radial, 
twist and mixed ones with different energies, and we observe whether the 
polarons are reflected, refracted or trapped. 

• Soft-inhomogeneity 

In this case, radial polarons moving along a chain of C-G base pairs 
are trapped by the soft A-T base pair as shown in Fig. El Trapping 
is caused by resonances between the polaron and the stationary state 
centered at the inhomogeneity. 
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Figure 4: Soft inhomogeneity and twist polarons. Transmission phenomenon 
due to the interaction between a moving polaron in a C-G chain with an A-T 
base pair 

However, when the polarons are activated by twist modes, the moving 
polarons are always transmitted. The inhomogeneity acts as a potential 
well, as shown in Fig.|31 Note that the polaron adopts briefly the shape 
of the ground state while passing through the inhomogeneity. 

• Hard-inhomog entity 

In this radial polaron moving along a chain of A-T base pairs in- 

teracts with the C-G inhomogeneity and is always reflected, as is shown 
in Fig. El This phenomenon is due to the impossibility of resonances, 
because the profiles of the radial components of the polaron and the 
stationary one located at the inhomogeneity are different. 

If the moving polaron is of the twist type, it is always transmitted as 
shown in Fig. IHl In this case the behavior is similar to the movement 
of a particle through a potential barrier, i.e, its velocity first decreases 
and then increases until its entering value. Note that again the polaron 
adopts briefly the shape of the ground state while passing through the 
inhomogeneity. 

In the mixed regime, the phenomena correspond exactly to the type of 
polaron activated. 
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Figure 5: Hard inhomogeneity and radial polarons. Reflection phenomenon 
due to the interaction between a moving polaron in an A-T chain with a C-G 
base pair 




Figure 6: Hard inhomogeneity and twist polarons. Transmission phe- 
nomenon due to the interaction between a moving polaron in an A-T chain 
with a C-G base pair 
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Finally, we have performed all the previous simulations introducing a 
small parametric disorder, in all cases the results are qualitatively similar. 
These can be summarized in the following table: 





Radial mov. 
regime 


Twist mov. 
regime 


Mixed regime. 
Radial pert. 


Mixed regime. 
Angular pert. 


A-T chain 


Reflection 


Trans. 


Reflection 


Trans. 


C-G chain 


Trapping 


Trans. 


Trapping 


Trans. 



The results are clear, the only possible candidate for charge transport 
mediated by polarons in an inhomogeneous chain are the twist modes. Ac- 
curate experiments are needed to find out good values of the parameters and 
for determining which kind of regime is possible in DNA and, therefore, if 
polarons may have a role in charge transport. 

Some preliminary numerical tests with a more realistic model where the 
ionization energy of a C-G base pair is lower than the A-T base pair by an 
amount of the order of 0.5 eV, show a trapping phenomenon, similar as the 
observed experimentally , where charge moves in an A-T chain by means 
of twist modes and reach a C-G base pair. In all other cases we have always 
observed a reflection phenomenon. This problem will be the object of further 
research. 

4 Conclusions 

We have considered a fully nonlinear, three-dimensional, semi-classical, tight- 
binding model for charge transport in DNA made out of identical base pairs 
except for a single one of different type, which we call a base pair inhomo- 
geneity. There are two types of this inhomogeneity, soft, composed of G-G 
base pairs with and A-T inhomogeneity and hard, with the complementary 
composition. 

In a previous work it has been described that in this system there 
exist two types of polarons, twist polarons and radial polarons, for which the 
electronic variables are coupled essentially to the radial or angular modes. 
The existence of these polarons depends both on the system parameters and 
on the form the movement of the polaron is activated. The properties of the 
two types of moving polarons are different. In general, the twist polarons 
are more robust with respect to parametric disorder, the polaron has lower 
velocity, and the activation energy is higher than in the radial movability 
regime. 

In an inhomogeneous chain with a single different base pair as a local 
inhomogeneity, we have observed that moving polarons activated by angular 
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perturbations are always transmitted by the inhomogeneity. If the polarons 
are activated by radial perturbations, we have never observed a transmission 
of the charge across the inhomogeneity. 

Some recent experiments on electron transfer in the DNA molecule 
show that electrons can migrate over long distances between a triplet C-G 
base pairs and a C-G base pair separated by a number n of A-T base pairs. 
Moreover, the triplet C-G base pairs acts as a sink for holes in the chain. 
In our model, decreasing slightly the ionization energy of the C-G base pair, 
as in real DNA occurs, we are able to reproduce the experimental results. 
Some more detailed numerical simulations are currently underway and will 
be published in due time. Also, these experimental techniques could be of 
application to contrast experimentally our results. 

In our model, we have considered a DNA chain in the vacuum. We have 
focused in the polaronic character of the charge carrier and its interaction 
with the chain. If the influence of the medium is not too strong, it is expected 
that the results would be similar as the ones obtained in the vacuum. In 
fact, they are in agreement with some experimental observations on DNA in 
aqueous solutions jHl]. The inclusion of thermal effects will be the subject 
of further studies, although in similar systems it has been shown that the 
main characteristics of the charge transport do not change very much when 
the temperature is considered [HK] . 

A Stationary polaron-like states 

In this section we expose the procedure followed for obtaining linearly stable, 
stationary localized states of our model given by Eqs. (jUHHI)- These solutions 
have been used in the previous section in order to generate mobile polarons. 
Since the adiabaticity parameter r is small enough, the fastest variables 
are the {c^}, with a characteristic frequency (the linear frequency of the 
uncoupled system) of order 1/r ~ 19, followed by the {r„} with frequency 
unity, and the {4>n} with ~ 0.11. We can suppose initially that r„ and 
(pn are constant, i.e., we use the Born-Oppenheimer approximation. For this 
purpose, we use a modification of the numerical method outlined in Refs. |2Hl 
EE]- We substitute in Eq. ^ Cn = exp(— i t/r), with time-independent 
$n's, and we obtain a nonlinear difference system = A^, with $ = 
($1, ...,$7v), from which a map $' = yl$/||y4$|| is constructed, ||.|| being the 
quadratic norm. 

Thus, using Eqs. (jBUHl), the stationary solutions must be attractors of the 



14 



map: 



r'n = -^|c„r (9) 



a 

K 

1, , 

aV 



— (c„c„_i + c„c„_i) + — (c„+ic„ + C„+iC„) ,(10) 

n+1 + 4>n-l) 

^/ ^ [(K + /c r^) c„ - (1 - g „) Cn+i - (1 - a rf^^^J c„_i] 

+ fcO c„ - (1 - a c„+i - (1 - a Cn-i}\\ ' 

where li' = d{r',(j)'). The starting point is a completely localized state given 
by c„ = 5ri,o, ■'"n = and 0„ = 0, Vn. Then the map is applied until 
convergence is achieved. In this way both stationary solutions and their 
energy E can be obtained. 

Firstly, we have analyzed the homogeneous and ordered case , i.e., = 
£"0, that arises in synthetic DNA (the constant Eq can be set to zero by means 
of a gauge transformation). As shown in Fig. [3 in a typical ground state the 
charge is fairly localized at few sites, and the amplitudes decay monotonically 
and exponentially with growing distance from the central site. The associated 
patterns of the static radial and relative angular displacements are similar 

We can introduce a certain degree of parametric disorder in the on-site 
electronic energy E^ by means of a random potential E^ G [—AE, AE] , with 
mean value zero and different interval sizes AE. In this case, the localized 
excitation patterns do not change qualitatively, as shown in Fig. |H| However, 
as the translational invariance is broken by the disorder, the localized excita- 
tion pattern is not symmetric with respect to a lattice site, which is different 
to the ordered case. Also, the localization is enhanced with the disorder, due 
to Anderson localization ^37ij . 

We consider now a chain with a local inhomogeneity due to the existence 
of a single base pair which is different to the other ones (without disorder). 
Our results show that the ground state centered in a site far from this local 
inhomogeneity is qualitatively similar to the obtained in the homogeneous 
chain case. Nevertheless, differences appear if we consider the stationary 
state centered in the inhomogeneity. In a C-G chain, with an A-T base pair 
as the local inhomogeneity, the shape of the ground state is qualitatively 
similar to the obtained in the homogeneous case, as shown in Fig. El In an 
A-T chain, with a C-G base pair as the local inhomogeneity, the static radial 
displacements are different, as shown in Fig. E3 As we have seen in the 
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Figure 7: Profiles of tlie ground state in tlie liomogeneous and ordered chain, 
(a) Wave function amplitudes |c„p. (b) Static radial displacements r„. (c) 
Static twists elongations 6*^71-1 
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Figure 8: Profiles of tlie ground state in tlie liomogcncous and disordered 
cliain. (a) Wave function amplitudes |C„|^. (b) Static radial displacements 
Tn- (c) Static twists elongations 9nn-i 
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radial displacements r„. (c) Static twists elongations 9nn-i 



previous section, this fact is determinant in relation with the transmission, 
reflection or trapping of a moving polaron by the inhomogeneity. 
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